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We consider the non-linear Poisson-Boltzmann theory for a single cylindrical or spherical macro- 
ion in symmetric 1:1, together with asymmetric 1:2 and 2:1 electrolytes. We focus on the regime 
where na, the ratio of the macro-ion radius a over the inverse Debye length in the bulk electrolyte, is 
large. Analyzing the structure of the analytical expansion emerging from a multiple scale analysis, 
we uncover a hidden structure for the electrostatic potential. This structure, which appears after a 
heuristic resummation, suggests a new and convenient expansion scheme that we present and work 
out in detail. We show that novel exact results can thereby be obtained, in particular pertaining to 
effective charge properties, in complete agreement with the direct numerical solution to the problem. 
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I. INTRODUCTION 

Colloidal suspensions are composed of large and often 
highly charged macromolecules (macro-ions or colloids), 
immersed in an electrolyte (“salt”) solution of mobile, 
positively or negatively charged micro-ions. These micro¬ 
ions move in a solvent that is in a first approximation 
regarded as a medium of uniform dielectric permittiv¬ 
ity. The system as a whole is assumed to be in thermal 
equilibrium at some inverse temperature j3 = l/ffceT). 

In the 1920s, Debye and Hiickel (DH) [I| proposed a 
linearized mean-field description of the bulk thermody¬ 
namics of Coulomb fluids which is adequate in the high- 
temperature region /3 —>■ 0. Some years earlier, Gouy 0 
and Chapman Q had established the nonlinear Poisson- 
Boltzmann (PB) mean-field treatment of the electric dou¬ 
ble layer, which served as a basis for the DVLO theory 
of colloidal interactions 0. 

A simple framework for studying the thermodynam¬ 
ics of colloidal suspensions at finite density is provided 
by the cell model • If however the concentration of 
colloids in the system is very low, in the first approxima¬ 
tion one can ignore their mutual interaction and consider 
the so-called infinite dilution limit. Colloids can be then 
studied as isolated mesoscopic objects of a given shape 
and bare charge, situated in a charged solution. This 
will be the viewpoint adopted here and for simplicity, we 
will address the homogeneous dielectric case with the di¬ 
electric permittivity of the colloid e' equal to that of the 
solvent e (no electrostatic image charges). 

The concept of (effective) charge renormalization, in¬ 
troduced by Alexander et al @ in the context of the 
PB cell model is simple to define in the infinite dilution 
limit: the renormalized/effective charge follows from the 
far-potential induced in the electrolyte [lol - [l^ . At large 
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distances from a charged body with bare charge Qbare, 
the electrostatic potential takes the same form as that 
obtained within the linearized DH theory, with a modi¬ 
fied prefactor Qeff which embodies the nonlinear effects 
of the PB theory, or of an approach which goes beyond 
the mean-field description. Within the nonlinear PB ap¬ 
proach, Qeff — Qbare for low values of IQbarel while Qeff 
saturates to a hnite constant when |Qbare| —>■ oo. In gen¬ 
eral, one expects that |Qeff| ^ IQbarel as a consequence of 
the nonlinear screening effect of the electric double layer 
around a colloid. For a monovalent 1:1 electrolyte the 
effective charge is indeed always smaller than the bare 
one. This is no longer true for asymmetric electrolytes 
which exhibit an overshooting effect [l^: there exists a 
(rather small) interval of Qbare where |Qeff| > IQbarel- 

Although the definition of an effective charge is unam¬ 
biguous within the nonlinear PB theory, it is not clear 
whether the far-potential behaves like the DH one in an 
exact description which goes beyond the mean-field as¬ 
sumption. The two-dimensional symmetric Coulomb gas 
of pointlike ig charges, interacting via the logarithmic 
Coulomb potential, is integrable in the whole interval of 
couplings 0 < P = < 2 where oppositely charged 

couples of charges do not collapse [l^. The concept of 
renormalized charge has been shown to be valid for a 
charged conductor wall [l^ , a pointlike guest charge M 
and for a guest charge with a small hard core [Tj, [^| 
which permits one to go beyond the stability thresh¬ 
old. It is interesting that for a guest charge with a 
small hard core at a finite temperature [l^, Qeff turns 
out to be a non-monotonous function of Qbare; the same 
phenomenon was observed also in Monte-Carlo [l^ and 
molecular-dynamics [1^ simulations of the salt-free cell 
model. Moreover, as Qbare goes to infinity, the effective 
charge does not saturate to a value but oscillates between 
two extreme (minimal and maximal) values. 

In this paper, we restrict ourselves to the definition of 
the effective charge within the nonlinear PB theory. For 
a charged infinite plane, one can obtain explicit results 
for the symmetric 1:1 and asymmetric 1:2 and 2:1 elec- 
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trolytes 0. For other asymmetric two-component elec¬ 
trolytes, the solution can be constructed implicitly, see a 
short review M- Realistic macro-ions are usually mod¬ 
eled as curved objects, namely cylinders or spheres of a 
given radius a. For such systems, the functional relation 
between Qeff and Qbare depends, besides the salt con¬ 
tent, also on the dimensionless parameter na where k is 
the inverse Debye (correlation) length of micro-particles. 
Two limiting cases are studied: 


• If KO <C 1, the analysis is difficult due to the 
counter-ion evaporation phenomenon, that may in 
the no salt limit be complete for spheres, or partial 
for cylinders, see and [l^. Yet, the cylindri¬ 
cal PB equation is Painleve integrable for the sym¬ 
metric 1:1 and asymmetric 2:1 and 1:2 electrolytes. 
This enables one to construct systematically the 
non-analytical Ko-expansion of the effective charge 


• If KO 1, the colloid radius is large compared to 
the Debye length k~^ and the plane geometry is 
a good reference for finding 1/{na) expansions of 
the effective charge. Shkel et al constructed an 
asymptotic large-distance expansion of cylindrical 
and spherical PB equations for the symmetric 1:1 
electrolyte. Using the method of multiple scales, 
they were able to derive the first two terms of the 
l/(Ka) expansion of the electric potential. Based 
on this work, analytical approximations were devel¬ 
oped in Ref. for the symmetric 1:1 electrolyte 
and in Ref. [1^ for the asymmetric 2:1 and 1:2 
electrolytes. 

The present paper concentrates on the large na limit, 
constructing expansions in l/(«:a) for the effective charge 
of cylindrical and spherical colloids. The m ultip le-scale 
method applied in the previous Refs. [l^[^,[^ is labo¬ 
rious and in practice does not allow to go to high orders 
of the l/(Ka) expansions. However, inspecting the struc¬ 
ture of these results suggests some re-summation can be 
performed, which in turn strongly points to a novel re- 
parametrization of the electric potential. Pushing further 
this idea, it appears that the algebra is conducive to a 
much easier derivation of the 1/(/ca) expansions. For the 
cylindrical geometry especially, the formulation ends up 
with a representation of each expansion order in terms of 
finite polynomials, which enables one to construct expan¬ 
sions to arbitrary high orders. As concerns the spherical 
geometry, we were able to go one order beyond the re¬ 
sults of [H, |2^; as a by product of the analysis, a 

divergence problem for higher-order terms indicates the 
change of the analytic behaviour of the expansion. For 
both cylindrical and spherical geometries and the three 
types of electrolytes (1:1, 1:2 and 2:1), the obtained an¬ 
alytical results for the coefficients of the l/(«:a) series 
expansions are tested against numerical resolutions. 

The article is organized as follows. In Sec. m we in¬ 
troduce basic formulae and definitions which are used 


throughout the paper. It is elementary and may be 
skipped by the reader familiar with the subject. The 
large-distance formalism of Shkel et al [2^ and the ensu¬ 
ing possible re-parametrization of the PB potential, are 
then explained in Sec. uni Our approach is presented for 
1:1 electrolytes in Sec. |IVl for 2:1 electrolytes in Sec. FVl 
and for 1:2 electrolytes in Sec. IVII A brief recapitulation 
and concluding remarks are given in Sec. I VIII 


II. BASIC FORMALISM 
A. Microscopic models 

In this paper, we shall consider the infinite dilution 
limit of colloids in suspensions, namely a unique colloid 
immersed in an infinite electrolyte solution. The system 
is formulated in the three-dimensional (3D) Euclidean 
space [ 13 . Three colloidal shapes are of interest: 

• The planar case, when the colloid occupies the half¬ 
space z < 0. The surface at z = 0 bears hxed sur¬ 
face charge density ere, where e denotes the elemen¬ 
tary charge and a has dimension [length]”^; with¬ 
out any loss of generality we assume that cr > 0. 

• The cylindrical geometry, when the colloid corre¬ 
sponds to an infinitely long cylinder of radius a, 
say along the z-axis, carrying a bare linear charge 
density Ae, A having dimension [length]”^. The 
system has a polar symmetry in the {x, y) plane. 

• The spherical geometry where the colloid is a 
sphere of radius a with center localized at the origin 
0 and carrying a bare charge Ze. Z is dimensionless 
and the system is radially symmetric. 

The space outside the macro-ion is filled by an infi¬ 
nite electrolyte solution. In general, it consists of M 
(= 2,3,...) types of mobile (pointlike) micro-ions v = 
1,2,..., M with positive or negative charges {qi/e}, where 
\q^\ is the valence of i^-species. The charged particles 
are immersed in a solvent which is a medium of uni¬ 
form dielectric permittivity e (in Gauss units, e ~ 80 
for water). They interact with each other and with the 
charged colloid surface via the standard Coulomb po¬ 
tential n(r) = l/(e|r|), which is the solution of the 3D 
Poisson equation 


A'(;(r) =-^(r). (2.1) 

e 

The system is in thermal equilibrium at the inverse 
temperature /3 = l/{kT); we denote by (• • •) the statisti¬ 
cal averaging over a thermodynamic ensemble. It is use¬ 
ful to introduce the so-called Bjerrum length /b = /3e^/e, 
i.e. the distance at which two unit charges interact 
with thermal energy kT. The macroscopic (i.e. ther¬ 
modynamically averaged over all possible particle con¬ 
figurations) density of (/-species at point r is defined by 
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n^{r) = {n^{r)}, where hi.(r) = * nu¬ 

merates the particles of species at spatial positions , <5 
denotes Dirac/Kronecker delta function/symbol and the 
hat in n means “microscopic” (i.e. for a given particle 
configuration). The charge density at point r is given by 

p(r) = (p(r)), p(r) = '^h^{r)q^e. (2.2) 


where a = 0 for the planar case, a = 1 for the cylindrical 
geometry and a = 2 for the spherical geometry. The cor¬ 
responding second-order differential equation (12.611 has to 
be supplemented by two boundary conditions (BCs), one 
at the contact with the colloid and the regularity one at 
an infinite distance from the colloid. The best way to de¬ 
rive these BCs is to consider the overall electroneutrality 
of the system. 


At large distances from an isolated colloid, i.e. in the 
bulk, the species number densities become homogeneous, 
n^{r) = rij/, and the requirement of the bulk electroneu¬ 
trality is equivalent to 

p{r) = p = ^n^q^e = 0. (2.3) 


B. Poisson-Boltzmann equation 


For a given charge density profile /o(rOi the mean 
electrostatic potential ip at point r is given by '0(r) = 
J dr'u(|r — r'|)/9(r'). The potential satisfies a counter¬ 
part of the Poisson equation dm, 

AV'(r) =-^p(r). (2.4) 

In the microscopic picture, the energy of an z/-particle 
at point r can be expressed in terms of the microscopic 
potential '!/'(r) as q^eipir) and the probability of finding 
the particle at r is proportional to the Boltzmann factor 
exp[—/3(3',ye'!/'(r)]. In a mean-field approach, one adopts 
this microscopic relation to the corresponding macro¬ 
scopic values, nj/(r) = rij/exp[—/3gi,e^(r)]; the normal¬ 
ization by the bulk value is consistent with the assump¬ 
tion that ■i/'(r) and its derivatives vanish in the bulk. This 
relation is exact in the high-temperature limit and only 
approximative for a finite temperature [2^. Considering 
it in the Poisson Eq. (1^ . one obtains a self-consistent 
PB equation for the mean electrostatic potential: 

AV^(r) = - ^ ^ n^q^ exp[-,5g,,e^/i(r)]. (2.5) 


In terms of the reduced potential (j){v) = Peilj{r) and 
the inverse Debye length k = the PB 

equation can be rewritten as 

A0(r) = - (2.6) 

All studied geometries are effectively one-dimensional 
problems. Let r be the distance from the plane in the 
planar case, the distance \/+ y'^ from the cylinder axis 
z or the distance \/+ y'^ + z'^ from the sphere center. 
The Laplacian for such symmetric case can be written as 

A t d 
r“ dr 


vj. . d^ ad 

t: = t::? + it:’ (2-7) 


dr 


dj-2 j. 


• The planar case: Integrating the ID Poisson 
equation 


d^ip{r) 

dr2 



( 2 . 8 ) 


over r from 0 to oo, we get 


47r 

Ip'(oo) — iIj'(0) = -/ drp(r). (2.9) 

e Jo 


The requirement of the overall electroneutrality 


pOO 

ae+ dr p{r) = 0 (2-10) 

Jo 


then implies the couple of BCs for the reduced po¬ 
tential 


4> {d) = —dnl'B.cr, lim cp'{r) = 0. (2-II) 

r—^oo 


• Cylindrical geometry: For a given charge den¬ 
sity profile of particles at r > a, the electroneutral¬ 
ity condition reads as 


poo 

Xe+ dr 27rrp(r) = 0. 

J a 


( 2 . 12 ) 


Multiplying the 2D Poisson equation 


1 d / d'0 

r dr \ dr 



(2.13) 


by r and integrating over r from a to oo, the con¬ 
dition of overall electroneutrality is consistent with 
two BCs for the reduced potential 

a(p'{a) = —2ZbA, lim r(p'(r) = 0. (2-14) 

r—)-oo 


• Spherical geometry: The electroneutrality con¬ 
dition reads 


/ 


Ze+ dr 47rr^p(r) = 0. 


(2.15) 


With regard to the 3D Poisson equation 


lA 

r^ dr \ dr / 



(2.16) 


the electroneutrality condition is equivalent to two 
BCs 


a?(p'{a) = —IbZ, lim r'^(p'{r) = 0. (2-17) 
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C. Effective charge 


with the BCs (12.111) takes the form 


The nonlinear PB equation (12.61) can be linearized by 
applying the expansion exp[—gj^(/)(r)] ^ 1 — q,^(j)(r). With 
regard to the bulk electroneutrality condition (j2.3l) . we 
arrive at the linear DH equation 

A0DH(r) = (2-18) 

whose form does not depend on the particular composi¬ 
tion of the electrolyte. This equation, supplemented by 
the appropriate boundary conditions, is solvable explic¬ 
itly for all considered geometries. 

The linearization of the potential Boltzmann factor is 
not adequate mainly in the neighbourhood of the col¬ 
loid, where the potential is large. On the other hand, at 
asymptotically large distances from the colloid the po¬ 
tential is infinitesimally small and the linearization pro¬ 
cedure is exact. We can say that the asymptotic PB 
solution satisfies the linear equation 

A4>(r) = |r| —>• oo. (2-19) 


0dh(7’) 


K 


(2.24) 


Let us choose Adh(o') = Airl-ga/K and f{Kr) = 
6 “”’’. Anticipating that the nonlinear PB poten¬ 
tial behaves asymptotically as 


cl){r) - A(CT)e-''^ (2.25) 

r—^co 


using the prescription (12.211) the effective surface 
charge density aese depends on the bare one ae as 
follows 


47r/B 

K 


CTeff = A(a). 


(2.26) 


• Cylindrical geometry: The linearized DH equa¬ 
tion 

0dh(^) “I—^dh(^) = «^^0dh(?’) (2.27) 

r 


Comparing with the DH equation (12.181) we see that the 
asymptotic forms of the PB and DH solutions are equiv¬ 
alent, up to position-independent prefactors: 

(('dh(?^) ^ ADH{Q)fiKr), 0(r) -- A{Q)f{Kr). 

r—^oo r—¥oo 

( 2 . 20 ) 

Here, the dependence of the A-prefactors on the ther¬ 
modynamic parameters of the electrolyte like na will not 
be explicitly indicated and Q is the bare charge charac¬ 
teristics of the colloid (the surface charge density ae for 
the plane case a = 0, the line charge density Ae for the 
cylinder a = 1 and the charge Ze for the sphere a = 2). 

The effective charge Qeff is defined as a function of the 
bare charge Q via the formula 

A{Q) = ADH((5eff)- (2-21) 


with the BCs (12.141) provides the solution 

<?i'DH(?') = — Ko(Kr), r>a, (2.28) 
KaKi[Ka) 

where Kq and Ki are the modified Bessel functions. 
Since KQ{Kr) ^ y/n({2Kr)e~'^'^ for asymptotically 
large r we can choose 


Adh(A) 


KaKi^na) 


XIb, 


f{Kr) 



(2.29) 

The prefactor to the asymptotic behaviour of the 
full PB potential 


(j){r) 


A(A) 



(2.30) 


In other words, QesiQ) is the effective charge in the lin¬ 
ear DH theory which reproduces the correct PB asymp¬ 
totic behavior; Qeff accounts for nonlinear effects, most 
prevalent close to the surface of the colloid. The nonlin¬ 
ear effects are negligible in the limit Q —0 and one ex¬ 
pects that QesiQ) '^Q-^o Q- In the opposite limit Q ^ oo 
one anticipates that the saturation value of the effective 
charge 


determines the dependence of the effective line 
charge density Agge on the bare one Ae as follows 

Aeff^B = _ KaKi{Ka)A{X). (2-31) 

V 27r 

• Spherical geometry: The linearized DH equa¬ 
tion 


QT^ = lim Qeff(Q) (2.22) 

Q-^oo 

is finite. 

Let us now assume that we know the prefactor A[Q) 
and derive an explicit formula for the effective charge for 
each of the three geometries. 

• The planar case: The solution of the DH equa¬ 
tion 

(2.23) 


—^dh(^) — ^^0dh(?’) (2.32) 

r 

with the BCs (12.171) has the solution 

r>a. (2.33) 

1 -I- Ata Atr 

It is natural to choose 

Adh(^) = finr) = —e” 

1 -|- Ata nr 


<^dh(^) = kVdh(?') 


(2.34) 
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Taking into account the expected asymptotic be¬ 
haviour of the PB potential 

(/)(r) - A{Z)—e-^^, (2.35) 

r—^oo tvT 

the formula for the effective charge as the 

function of the bare charge Z reads as 

(2.36) 

a na 

In other words, Z^s is the value that should be 
plugged in the right-hand side of (12.331) . so that 
the latter formula provides the correct far-held of 
the non-linear solution to the PB equation. By 
construction thus, ZeS = Z when the PB theory is 
linearizable, which is the case for Z —0. 


D. Explicit results for the planar case 


• Asymmetric 2:1 electrolyte. In this case, the 
positively charged colons to the surface have gi = 2 
and the negatively charged counterions have (72 = 
— I. The bulk electroneutrality requires that rii = 
n/3 and 712 = where n is the total par¬ 

ticle number density. The inverse Debye length 
K = -v/Stt/bu. The PB equation 


(j)''{r) = 


3 1 


e0(r) _ g-20(r) 


has been solved by Gouy 0: 


(2.42) 


(f>{r) = In 


1 + 


A(cr)e-'"’’ 



(2.43) 


The potential behaves at large kt as (12.251) . The re¬ 
lation between A and a follows from the BC (12.111) 
at r = 0: 


The planar case is solvable explicitly only for specific 
types of two-component electrolytes. 

• Symmetric 1:1 electrolyte. We have two types 
of particles with (reduced) positive qi = 1 and neg¬ 
ative 92 = —1 charges. Denoting by n the total par¬ 
ticle number density, the requirement of the bulk 
electroneutrality is equivalent to ni = 77,2 = n/2. 
The inverse Debye length is given by k = \/47rlBn. 
The corresponding PB equation 


47r/Bq- ^ 36A(6 -f A) 

K (6-A)(A2-h24A-k36)' ^ ^ 

From the three A-solutions of this cubic equation 
we take the one which goes to zero in the limit 
cr —>■ 0. In the saturation limit a —>■ 00 , we have 
Asat = 6 and therefore 

= 6. (2.45) 


(f>" (r) = sinh 
has the explicit solution 


(j){r) = 2 In 


1-f 


AM, 


1 _ AMe 

r 4 e 


— Kr 


— Kr 


(2.37) 


(2.38) 


• Asymmetric 1:2 electrolyte. Now the coions 
have <71 = 1 and the bulk number density ni = 
2nj3, while the counterions have <72 = ~2 and 77.2 = 
77/3. As before, k = ^/SttIbu. The PB equation 


cj)"{r) = 




(2.46) 


which indeed behaves at large kt as predicted by 
formula (12.251) . The relation between A and a is 
yielded by the BC (12.111) at r = 0: 


A{a) = 


SttZbct 


1 


i + Vi + (^)^ 


(2.39) 


In accordance with the relation (12.261) , the effective 
charge density is given by 

2a 

tTeff = I ■ (2.40) 

l + Vl + (^) 


It has the correct behaviour a^s ^ cr in the limit 
tj —^ 0. In the saturation cr —>■ 00 limit, we have 


has the explicit solution 


4>{r) = — In 


1 - 


A(cr)e 


^1 -t 


(2.47) 


The relation between A and ct, following from the 
BC (12.111) at r = 0, takes the form 




36A(6 - A) 


(6-f A)(A2-24A-k36)’ 


(2.48) 


The physical A-solution goes to zero in the limit 
cr —^ 0. In the saturation limit cr —>■ 00 , Asat =6(2— 
•\/3) is the smaller root of the quadratic equation 
A2 — 24A -I- 36 = 0 and we arrive at 


47r^B 

K, 


—sat 


= 4. 


(2.41) 


=6(2-73). (2.49) 
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III. AN ASYMPTOTIC EXPANSION FOR 1:1 
ELECTROLYTE 


For the symmetric 1:1 electrolyte of total particle num¬ 
ber density n and k = y/AirlBn, the general PB equation 
for the electrostatic potential (j) reads as 

-|- —(j)'(x) = sinh^(x), (3.1) 

X 


where x = nr is the reduced distance, a = 1 for the 
cylindrical colloid and a = 2 for the spherical colloid. In 
the paper [ 2 ^, an asymptotic large-x expansion of the 
potential has been obtained in the following form 


(j){x) = 


3,0/2 

a,-Sx 


X' 


,3a/2 


IQO 


A 


Ao2 

~Ar 


10 


^01 

X x^ 

^11 ^12 

' ' ^ 

X 




+ 


g-(2j + l)T 

^(2j+l)a/2 


A 30 + ^ 


Aj2 



Here, Aqo = H is the crucial prefactor to the lead¬ 
ing large-distance asymptotic. The other prefactors Ajk 
scale with A like Aj^ = where the coefficients 

(numbers) {djk} fulfill a rather complicated recursion 
which enables one to generate systematically the coef¬ 
ficients {ajk} with small values of the {j,k) indices [25j |. 

The first j = 0 row of the asymptotic formula (13.2p . 
which consists in the exponential e““ multiplied by an 
infinite inverse-power-law series, corresponds to the lin¬ 
ear DH approximation. The next rows are exponentially 
smaller and smaller corrections to the DH result. 

For our purposes, it is more important to concentrate 
on columns. Let us introduce the new variable 


A 

= (3-3) 

where A (and thus t) implicitly depends on the bare 
charge (A or Z depending on the geometry), and na; we 
shall usually omit in the notation these functional de¬ 
pendences or write only the relevant ones. In the BC, 
the first derivative of the potential is taken just at the 
surface of the colloid, i.e. at x = na. The value of t is, 
in general, not small at the colloid surface, even in the 
limit of interest na ^ 00 . Note that the given column 
differs from the previous one basically by the factor I/x. 
The corresponding surface factor l/{Ka) is small in the 
limit Ko —>■ 00 . This permits one to generate a systematic 
expansion of the prefactor A in 1 / (ku) powers by taking 
successively column by column. 

The asymptotic expansion (13.21) can be rewritten in the 
variables t and I /x as follows 


j 

0 

II 

k=i 

k =2 

CO 

II 

II 

0 

4 

1 

~2 

9 

25 

75 

28 

3675 

213 

1 

4 

3 

71 

3215 

79521 

3 

2 

25 

3 X 28 

213 

9 

4 

3 

89 

1485 

114071 


5 

2 

25 

28 

213 

Q 

4 

3 

107 

1965 

157237 


7 

2 

25 

28 

213 

4 

4 

3 

125 

2517 

211275 


9 

2 

25 

28 

213 


TABLE I. Cylindrical geometry. The coefficients Bjk {j, k = 
0,1, 2, 3, 4) appearing i n (13.51) for the large-distance expansion 
(EH) obtained in Ref. 


are the numbers which will be explicitly available for 
small values of (j, k) indices. 

We shall also need the asymptotic large-x expansion of 
the modihed Bessel functions 


K^{x) 



1 ^(5 1' k) 

(2x)'=fc!F(i-b^-A:) 


(3.6) 


for z/ = 0 , 1 . 


A. Cylindrical geometry 

For the a = 1 cylindrical geometry, we have the pa¬ 
rameter 


t(x) 


A(A)e-^ 
4 ^/x 


(3.7) 


The coefficients Bjk with j,k = 0,1, 2, 3,4 are summa¬ 
rized in Table HI The first row of the table corresponds 
to the DH result p.28|l . Indeed, writing in (12.281) the 
large-x expansion of the modified Bessel function (13.61) 
with 12 = 0, the corresponding coefficients are found to 
be 


00 ^ 

j,k=0 

where the coefficients 

_ 

j ^jk 


(3.4) 


Bok = 4- 




Ok 

~A~ 


= 2 


2-k ^(5 + k) 


k\r{A-k)' 


(3.8) 


the first few of which read 4, —1/2, 9/2® etc. 

In accordance with our strategy, let us consider the 
first column of Table HI From the first few coefficients 
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BjQ we can “guess” 


Bjo - 


4 

2j + l 


(3.9) 


and suggest that this formula holds for all j = 0,1 ,.... 
If this is true, the potential (13.41) is given, in the leading 
Xjx^ order, by 


OO 


<))(x) =4^ 

f=o 


^ 2 ^ + 1 
27+1 



(3.10) 


This result is identical to the planar one (12.381) under 
the only proviso that the cylindrical t (IS21) differs from 
the corresponding planar function Ae~^ 14 by the factor 

1/^/i. 

We can go further and sum the contributions of the 
second column of Table U to determine the potential up 
to the 1/x order. We anticipate that 


i?oi - 




(3.11) 


The potential (13.41) is then given by 
1 t 


4i{x) = 2 In 


= 2 In 


1 - t 
1 t 
1 - t 


1 3 

- 1 -^ 

2x 2x 1 — 

t(l + 2^2) / 1 


4x 


l-t 1+t 


(3.12) 


The l/x correction consists of 1/(1 — t) and 1/(1 + t) 
terms. They arise naturally from a “renormalization 
ansatz” 


with 

/(x) =t(a;)|l-^[l + 2t^(a;)]| . (3.14) 


It is easy to verify that this ansatz coincides with the 
original equation (13.121) to the order 1/x. 

The guessing of the coefficients becomes harder when 
considering the next columns. In the following section, 
we shall show how to generate systematically the whole 
infinite series of coefficients in a straightforward way. At 
this stage, we restrict ourselves to the preliminary result 

(Eni), (EH. 

For the solution of type (13.131) . the BC (12.141) at r = a 
can be expressed as 

This relation determines the function A(A). In the sat¬ 
uration limit A —>■ OO we have /sat(Ko) = 1, as indicated 
by Eq. (13.151) . Consequently, 


j 

k=0 

k=l 

k=2 

k=3 

k=4 

0 

4 

0 

0 

0 

0 


4 


7 

15 

155 


3 


2 

2 

23 


4 


14 

35 

2327 


5 


3 

3 

32 X 23 

Q 

4 

9 

205 

863 

40621 


7 


32 X 22 

33 X 2 

33 X 25 

4 

4 

2 

2413 

90271 

11313029 
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5 X 32 X 23 

5 X 33 X 25 

52 X 33 X 28 


TABLE II. Spherical geometry. The coefficients Bjk {j, k = 
0,1, 2, 3, 4) appearing i n (13.51) for the large-distance expansion 
(1^ obtained in Ref. [^. 


Performing the large-KO expansion in this formula, we 
obtain 


tsatina) — 1 + -- [-••• = 

SKa 4 \/Ka 

where Agat = A(A —>■ oo). Using that 


(3.17) 


^iM = 1 + ^ 


2Ka 




(3.18) 


the formula for the effective charge ()2.31|1 implies that its 
saturation value has the large-Aca expansion of the form 


\ S3itl 1 _ 

^eff — 


1 




KaKi{na)Asa,t 


= 2Ka - 


(3.19) 


This result agrees with the previous finding of Ref. 


B. Spherical geometry 

For the a = 2 spherical geometry, we have 

i(x) = j —. (3.20) 

4 X 

The coefficients Bjk with j,k = 0,1, 2, 3,4 are summa¬ 
rized in Table IE 

The first row of the table corresponds to the DH result 
(12.281) . namely i?ofc = 4(5ofe. The first column of Table HD 
is identical to the first one of Table HI so that 


1 - [1 ■ 

Sna 


2 tL(Ka)] 


Bjo - 


4 

2j + l 


^sat (^n) 


= 1. (3.16) 


(3.21) 


























































and, in the leading 1/x^ order, the potential is given by 

1 + 1 ' 


4>{x) = 2 In 


1 - t 


(3.22) 


Again, this result is identical to the planar one ()2.38ll if 
the planar t is multiplied by 1/x. 

It is easy to guess the second column of Table [iTl 


Bqi = 0, Boj = -2 for j > 1. 
The potential (13.41) is then given by 

'l + t\ ^ 1 1 


= 2 In 


1 - t 


1-t 1+t 


(3.23) 


(3.24) 


The potential is transformable to the renormalized form 
of type (|3.13|) with 


fix) = tix) 




(3.25) 


IV. GENERAL FORMALISM FOR 1:1 
ELECTROLYTE 

The above section motivates us to search for the poten¬ 
tial 4>ix) in the ansatz form (13.131) which is nothing but 
the redehnition of the potential in terms of a new func¬ 
tion fix) with simpler expansion property, as we will see 
later. Inserting the ansatz (13.131) into the PB equation 
we obtain after some simple algebra the following 
differential equation for the /-function: 


fix) + -fix) - fix)\ [1 - fix)] 

-2f{x) [fix) + fix)] [fix) - fix)] = 0. (4.1) 


With regard to the results of the preceding section, we 
expect that the /-function can be written as 


fix) = tgix), 


(4.2) 


where the x-dependent function t is defined in (13.3p and 
the g-function is expected to have the following large-x 
expansion 


As before, the saturation limit Z —>■ oo corresponds to 
/sat(Ka) = 1- Consequently, 


4at(Ka) 


1 - 


2 /ta 




= 1 . 


(3.26) 


The large-KO expansion of tgatina) then reads 


tsatfa) = 1 - 1 - -— 
ZKa 


Asat e-- 
4 na ’ 


(3.27) 


where Agat = A(Z —>■ oo). The formula for the effective 
charge (|2.36p then leads to 


vsat^B _ l + Ka 

i 


a KQ 
= Ana - 1 - 6 - 1 - 


(3.28) 


which coincides with the finding of Ref. [1^. It is in¬ 
structive to compute the corresponding effective surface 
charge (7®^ = /iAnaf and likewise for the cylinder, 

CTgg = Agff/(27ra). In doing so, it appears that Eq. (13.281) 
and Eq. (13.191) bear the same information. Indeed, intro¬ 
ducing the curvature C = 1/a for cylinders and C = 2/a 
for spheres, both can be written, when phrased in terms 
of the effective surface charge, as 


_sat 

^eff 


/-(i + |£) 

ttIb V 


(3.29) 


To dominant order, that is when C = 0, we recover as 
expected the planar result of Eq. (12.411) . This suggests 
that to dominant plus sub-dominant order, the effective 
charge only depends on the curvature, irrespective of fur¬ 
ther geometrical details [1^ . We shall see that this “uni¬ 
versality” is broken by higher order terms in C^. 


OO ^ 

gix) = l + ^-fe<?fc(0- 


(4.3) 


k^l 


In the DH limit t —>■ 0, / is going to 0 as well and from 
f|3.13p we can write that (j){x) ^ 4/(x). Using the original 
expansion (ixii) . we identify 

9ki0) = ^ for all fc= 1,2,.... (4.4) 

From the definition of t in Eq. (1331) . we have that 



(4.5) 


Inserting the representation (14.2p into Eq. 63) and us¬ 
ing these relations, we obtain the differential equation for 
the g-function: 


fix) - 2g'ix) + 


[l-fg^ix)] 


2fgix) g'ix) - ^df) gf) - gfx) 


= 0 . 
(4.6) 


We emphasize that a prime here means the total deriva¬ 
tive with respect to x, including the function t(x). Like 
for instance, within the representation (14.31) we have 


g'ix) = -^^ \kgkit) + |t<?;(t)] -^ -ktg'kif 






(4.7) 

and so on. The point is that the total derivative with 
respect to x keeps the power-law expansion in 1/x where 
each term is multiplied by a function depending on t 
only. This two-scale method permits us to determine 
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recursively the functions {gkit)}- Setting to zero the 
coefficient to 1 /x, we obtain a differential equation for 
gi{t), supplemented by the BC gi(0) = Bqi/A deduced 
from (gH). Setting to zero the coefficient to l/x^, we 
obtain a differential equation for 32 (t) which involves 
also the known function qi(t), supplemented by the BC 
52 ( 0 ) = -B 02 / 4 , etc. 


A. Cylindrical geometry 


We start with the cylindrical geometry for which a = 
1. The BC (14.4p implies 


fffe(O) = 


1 r(i + fc) 

2^= k\T{^ - k) 


(4.8) 


for all fc = 1, 2,.... 

Setting to zero the coefficient to 1/x in (14.61) . gi{t) 
must obey the differential equation 


t(l — t^)gi{t) + (3 + t^)gi{t) + 2t — 0. ( 4 - 9 ) 


etc. In general, gkit) with k = 1,2,... turns out to be 
a finite polynomial of the fcth order in the variable 
This special and convenient property of the ^fe-functions 
is present exclusively for the case of the cylindrical ge¬ 
ometry and the symmetric 1 : 1 electrolyte. 


-0.38 r 


-0.4 - 

$ 

"7 -0.42 - “ ® 
a 

^ ■ -0.44 - 

i 0-75 

~ -0.46 
a 


■'vD 


0 


-0.48 - 


-O.Si; 






0.02 


0.04 


0.06 


0.08 


0.1 


(Kfl)' 


The general solution of this equation, obtained by using 
the ‘variation of constants’ method, reads as 


giit) = Cl + C 2 t^ 


(1 -b 4 c 2) (Int -b 


1 

4^2 


(4.10) 


The integration constant C 2 is determined by the regu¬ 
larity of giit) as t —>■ 0 as follows C 2 = —1/4. The BC 
5i(0) = —1/8 leads to ci = —1/8. The consequent 

giit) = -\ii + 2e) (4.11) 


is in full agreement with the previous result (13.141) . 

To find 52 ( 0 , we set to zero the coefficient to 1/x^ in 
(14.61) which, together with the knowledge of gi(t), implies 

t(l-t 2 ) 5 ''(t) + (3+t2)5'(t)-pi(2t4-6t2-i5) = 0 . (4.12) 

The general solution of this differential equation is 

»(*) = ”■ + + T 6 (f - '■"») + i) ■ 

(4.13) 

The regularity of g 2 {t) at t = 0 implies C 2 = 15/32 and 
the BC 32 ( 0 ) = 9/128 fixes ci = 9/128. Thus we arrive 
at 

, , 32 3 X 5 2 1 4 

92{t) - ^ + • (4-14) 


In the same way, e.g. by using the symbolic language 
Mathematica, we get 


53 (i) 

giit) 


3 X 52 
210 


163 X 3 2 _ ^ ,4_^ ,6 

29 27 43 ’ 


(4.15) 


3 X (5 X 7)2 
215 


3 X 53 X 59 

212 


r -b 


34 X 17 

211 





1 8 
J_ 

+ 44 ^ 


(4.16) 


FIG. 1. Asymptotic expansion plots for cylinders in a 1:1 
salt. The main graph shows Ka(Agg Zs — 2/ta — 3/2) as a func¬ 
tion of l/(Ka). The dashed linear curve is for the first order 
correction predicted by Eq. (14.181) . — 1/2 -b 27/(32Ka), while 
the parabolic dot-dashed curve is for —1/2 -b 27{Ka)~^/32 — 
131(Ka)“3/64, The symbols correspond to the saturated ef¬ 
fective charge obtained from the numerical solution of Eq. 
(ED with Q = 1, and have been computed for na = 10, 
30, 50 and 100, as can be read on the (c-axis. Inset : 
(Ka)3[A|gZs — 2Ka — 3/2 — 0.5/(Ka)] versus l/(fi:a). The lin¬ 
ear dashed curve displays the prediction of Eq. (ED for the 
quantity of interest, namely 27/32 — 131(/ta)“i/64. 


Having fix) (14.21) with g(x) (14.31) truncated say at 
k = A, the crucial function A(A) is again determined 
by the relation (13.151) . In the saturation limit A —>■ 00, 
the requirement /sat(7«a) = 1 implies the iteratively gen¬ 
erated large-Ko expansion 


, , 3 1 

tsat(Ka) = 1 + - 

2'7 na 


5x7 1 

27 (Ka)2 


7369 X 5 1 

215 (Ka)4 + ■ ■ ■ 


53 X 32 1 

213 (Ka)3 

(4.17) 


which goes beyond the previous one (13.171) . Using the 
connection (13.191) . the corresponding saturation value of 
the effective charge exhibits the large -ko expansion of the 
form 


Ks^b — 2Ka + - — 


1 1 
2 na 


131 1 

64 (Ka)3 


+ 0 


27 1 

32 (Ka)2 



(4.18) 


Notice that in spite of the complicated large-Ata expansion 
of tsat(6(a), the corresponding expression for A|g is rel¬ 
atively simple. The numerical checks of the coefficients 
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to the l/(Ka), l/(Ka)^, ... terms are presented in Fig. 
HJ Such a comparison poses the difficulty that the ef¬ 
fective charge Agg be known with high precision in the 
limit where the bare charge A diverges. To this end, the 
Poisson-Boltzmann equation m is solved numerically 
for a series of increasing bare charges, at a given value 
of Ka. The effective charge Agg is extracted from the 
far-held behaviour of the potential, that reads 


The regularity of gi{t) as t —>■ 0 hxes C 2 = —1/2 and the 
BC 5 i( 0 ) = 0 implies Ci = 0. Thus we have 

g,{t) = -^t^ (4.23) 

which agrees with the result (13.251) . 

Setting to zero the coefficient to Xjx^ in (14.61) . we ob¬ 
tain the differential equation for 32 (t): 


r ^ 00 . (4.19) 

KaKi(Ka) 

A ‘hnite-charge’ scaling analysis is subsequently per¬ 
formed: it is indeed straightforward to show that Agg — 
A|g vanishes as 1/A, when A —>■ 00 . In practice, the above 
linear regime in 1/A is well reached whenever A > 10®. 
The saturated values thereby obtained are shown by the 
symbols in all the graphs displayed. Once A|g is known, 
inspecting its behaviour as a function of Ka, as performed 
in Fig. [TJ allows form a stringent test of the analytical 
prediction. The main graph reveals that beyond the dom¬ 
inant behaviour in 2Ka + 3/2, the next correction to A^^ 
is —l/( 2 Ka), since the data shown extrapolate to — 1/2 in 
the limit l/(Ka) —>■ 0. Besides, the next term predicted 
with prefactor 27/32 brings significant improvement at 
large although finite Ka (see the linear dashed line in the 
main graph). Yet, some (negative) curvature can be in¬ 
ferred from the symbols shown and indeed, including the 
next term with prefactor —131/64 further enhances the 
agreement. The inset offers a direct proof that expression 
(14.181) . with all terms, is a very plausible expansion for 
the saturated effective charge. Note that the quadratic 
dashed-dotted curve in the main graph and the dashed 
line of the inset convey the same information, in a differ¬ 
ent visual setting. 

We would like to emphasize that our method of gen¬ 
erating the large-Ko expansion is technically very simple 
and we can generate within few seconds by using Math- 
ematica also the next higher-order terms of the series 

mB- 


B. Spherical geometry 


For the spherical a = 2 geometry, the BC (14.41) implies 


5fe(0)=0 (4.20) 


for all A: = 1, 2,.... 

Setting to zero the coefficient to 1/x in (14.61) . gi{t) 
fulfills the differential equation 


til-t^)g'^it) + {3 + t^)g'^{t)+ti2t*-7t^-7) = 0. (4.24) 

The requirement of regularity as t —>■ 0 and the BC 
52 ( 0 ) = 0 imply the solution 


771 

g^(t) = + - Y 


2 ^ n2(n2 - 1)' 


(4.25) 


The function 52 (t) is an infinite polynomial in with the 
convergence radius t < 1. In particular, 


..(!) = 2--. 


(4.26) 


0.375 


0.35 

a 

-5 

0.325 

y 

0.3 


0 0.01 0.02 0.03 0.04 0.05 

(Kfl)'^ 

FIG. 2. Asymptotic expansion plots for charged spheres in a 
1:1 salt (main graph), together with asymmetric 1:2 and 2:1 
salts (insets). The main graph shows «;a(Z||‘Zs/a — 4«:a — 6) 
as a function of l/(«:a). The dashed linear curve is for the line 
that passes through the two points associated to the largest 
values of Ka (50 and 100). The arrow is the predicted asymp¬ 
totic limit 7r^/3 — 3. The symbols correspond to the satu¬ 
rated effective charge obtained from the numerical solution 
of Eq. (13.11) with q = 2. The 2:1 inset shows the coun¬ 
terpart, Ka(Z|g‘lfl/a — Gna — 7) where the limit is predicted 
to be ~ 3.844. The 1:2 inset is constructed similarly from 
Eq. (16.2111 . for which the arrow indicating the intercept is at 
-0.089. 



t{l-f)g'{{t) + {3 + t^)g[{t)+it = 0. 


The general solution of this equation reads 


(4-21) In the saturation limit A —> c», the requirement 
fsati^a) = 1 implies the iteratively generated large-Ata 
expansion 


gi{t) = Cl + C 2 t^ 


2(l-h2c2) 



(4.22) 


, , 11 
tsat(Ka) = 1 + -- 

z na 


5 TT^ A 1 
4 ~ 12/ («;a)2 


(4.27) 
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Here, we used that 52 ( 1 ) is given by (14.261) . We cannot go 
beyond the indicated order because the next term needs 
the iteration with the diverging value of 52(1 + 1/(2ko)). 
This indicates that the next-order singular term has a 
form different from l/^na)^. 


V. 2:1 ELECTROLYTE 

For the asymmetric 2 : 1 electrolyte in contact with 
the cylindrical a = 1 or spherical a = 2 colloids, the PB 
equation in the reduced distance x = nr takes the form 


Using the prescription (I3.28L the corresponding satu¬ 
ration value of the effective charge exhibits the large-KO 
expansion of the form 


0"(x) + -ct>'{x) = i . 

X o L 


- 24 >{x) 


(5.1) 


Being motivated by the exact planar solution (I2.43L 
we search the electrostatic potential in an ansatz form 


/7sat (b 

a 


= Ana -1-6-1- 




(4.28) 


(l)(x) = In 


1 


_6f{x)_) 
[1-/(2^)]^ ■ 


Introducing 


t{x) 


A{X) e-^ 

6 ’ 


(5.2) 


(5.3) 


The numerical check of the prefactor to the third l/(/ca) 
term is presented in Fig. [5J Indirectly, the plot also 
assesses the correctness of the dominant terms 4^0 -(- 6. 

I 


f{x) ^ t{x) is the leading large-distance form. Inserting 
the ansatz (15.21) into the PB equation (15.11) . we obtain 
the differential equation for the /-function: 


- f{x) [1 - f{x)] [fix] + 4f{x) + 1 ] 

-2 [f{x) + 3 fix) + 3fix) - 1 ] [fix) -b fix)] [fix) - fix)] = 0. 


We again expect that the /-function is expressible as a series 


(5.4) 


fix) = tgix), 


00 ^ 

gix) = ^ + J27k9kit)- 




(5.5) 


The g-function satisfies the differential equation 


g"ix) - 2g'ix) + 


[1 - fg'^ix)] fg'^ix) + Atgix) + l] 


-2t fg^ix) + 3fg^ix) + itgf) - l] gfx) - ^gf) (f + gix) - gfx) 


= 0 . 


(5.6) 


A. Cylindrical geometry 


For the cylindrical a = 1 geometry, setting to zero 
the coefficient to 1/x in (15.61) . gi(t) obeys the differential 
equation 


t(l +4t — — tfg'fit) 

-|-(3 -|- 8t -|- 12t^ -|- f)g[it) 

+2i-l + 3t + 3f+t^) =0. (5.7) 


The solution of this equation with the BC gi(0) = —1/8, 
regular as t —>■ 0, is found to be 


gi{t) 





A 1 


Similarly, we obtain 


52 (i) 


gft) 


2543 11 5 2 1 3 1 4 

2^x32^2x32 25^ ^32 ^24x32 

^2x51 5x7 1 ^23 1 

^ 3 l + t“ 2x32(1+ t)2 +^(l + t)3> 

(5.9) 

5 X 41 X 659 , 72 ^ 617 ^2 17 ^3 

210x33 24 X 32^” 29 X 3^ ~ 22 x 33^ 

13 X 19 4 ^ 1 ,5 _ 1 .6 

27 X 33 23 X 32 26 X 33 

193 1 , 11173 1 

~ 2 X 32 1 -b t 23 X 33 (1 -bt)2 

4649 1 ^ 43 1 2® 1 

~ 24 X 33 (1 -bt)3 ¥ (1 -bf)4 “ ¥ (1 -bt)3’ 

(5.10) 


3 1 + t 


(5.8) 
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etc. In general, gkit) (fc = 1, 2,...) is the sum of two finite 
polynomial, one of the 2 A:th order in t and the other of 
the {2k — l)th order in 1/(1 + t). 

In the saturation limit A —>■ oo, the requirement 
/sat(Ka) = 1 implies the iteratively generated large-KO 
expansion 


tsat{Ka) = 1 + 


5 1 

2^ X 3 Ko 


5 X 41 1 

27 X 32 {Kay 


12583 1 _ Al,,t 

2^0 X 32 {na)^ 6 -y/KO 


(5.11) 


Considering the BC gi(0) = 0, the solution is 


9i{t) 


8 4 1 2 8 1 

- 1 --. 

336 31+t 


(5.14) 


Setting to zero the coefficient to 1/a:^ in (15.6L we ob¬ 
tain the differential equation for g 2 {t) of the form 


P{M{t) + Q{t)g!,{t) + R{t) = 0, (5.15) 


where 


Using the prescription (13.191) . for large values of na the 
saturation value of the effective charge behaves as 


A 


sat 

eff 




Zna -\ -h 

4 


5 1 

12 KO 


703 1 

{naf 


+ 0 


1 


{naY 

(5)12) 


The first two terms of this series have been obtained in 
Ref. M- The numerical checks of the coefficients to 
the l/(Ka) and l/(Ka)^ terms are presented in the main 
graph of Fig. [3l 
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FIG. 3. Asymptotic expansion plots for cylinders in 2:1 and 
1:2 salts. The main graph for the 2:1 case shows fi:a(Aeff is — 
Sna — 7/4) versus !/(««). The linear dashed curve displays 
the prediction embodied in Eq. (15.1211 . 5/12 — 703(Ka)“^/192. 
The inset displays expression (16.131) as a function of !/(«:«). 
The dashed linear curve is for the prediction of Eq. (16.1211 . 
—0.295 -I- 0.700/(K:a). As in previous figures, the symbols 
show the data obtained from solving numerically the Poisson- 
Boltzmann equation. 


P{t) = 9t{l -h tY{l + 3t- 3Y - Y), (5.16) 

Q{t) = 9(1 -h t)4(3 + 8t + UY + Y), (5.17) 

R{t) = 78 -h 51t -h 8AY + 252 Y - I 912 Y 
-2560t® - 1860t® - 430t^ -h 314t® 
-hl89t®-h32U°-h2t“. (5.18) 


The function 52 (i) is an infinite polynomial in t which 
diverges for alH > 1. The value of 52 at t = 1 is of our 
primary interest. Denoting 52 ( 1 ) = Al, by using Mathe- 
matica it can be shown that 


— i ln(2 — Vi) ln(3 — Vs) 

fj 


1 

~6 L 


ln(3 -h V^) 

12 ~ 5/3)^ — -Li 2(—2 + Vs) 


~ -0.27962, 


(5.19) 


where lji 2 {z) is the polylogarithm function defined by 


00 ^ 

= 1^1 ^ 1 - ( 5 . 20 ) 

The saturation condition fsat {kg) = 1 implies the 
large-Ko expansion 


isat(Ka) — 1 -|- —-h 

6 na 



1 

1^'" 


A „-Ka 
■^sat ^ 

6 Ka 


(5.21) 


With respect to the prescription p.28|) . the saturation 
value of the effective charge behaves for large values of 
Ka as follows 


B. Spherical geometry 

For the spherical a = 2 geometry, setting to zero the 
coefficient to 1/a; in (15.61) . gi{t) fulfills the equation 

t(l -I- 4t — — Y)g"{t) 

-l-(3 -|- 8 t -l- 12t^ -|- Y)g\{t) 

+A{-l + 3t + 3Y+ Y) = 0. (5.13) 


+ 7 + ± + o 

a \ 6 J Ka \Ka J 

~ 6 Ka-f 7-h3.844—-hof — V (5.22) 

Ka \Ka J 

The first two terms of this expansion are in full agreement 
with the result of Ref. dl. The prefactor to the l/(«:a) 
term is checked against the numerical resolution in Fig. 
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VI. 1:2 ELECTROLYTE 


With regard to the exact planar solution (j2.47|) . the elec¬ 
trostatic potential is searched in an ansatz form 


For the asymmetric 1 : 2 electrolyte, the PB equation 
takes the form 


(j){x) = — In 


_6f(x)_] 

[l + /(^)f J ■ 


Introducing 


t{x) 


A{X) e-^ 

6 x°‘/'^ ’ 


( 6 . 2 ) 


(6.3) 


'{x) + %'{x) = \ 
X 3 


q 24 ,{ x ) _ g-0(x) 


f{x) ^ t{x) in the leading large-distance order. 

Inserting (16.21) into the PB equation (16.IL the /- 
( 6 . 1 ) function obeys the differential equation 

I 


f"{x) -k -fix) - /(a;)J [1 - p{x)] [f{x) - 4/(a;) -H l] 

-2 [f{x) - 3 fix) + 3fix) + 1] [fix) 4- fix)] [fix) - fix)] = 0. 


Writing the /-function as a series 


(6.4) 


fix) = tgix), gix) = l + '^—gkit), 


(6.5) 


k^l 


the ^-function satisfies the differential equation 


g"ix) - 2g'ix) + 


[l - fg'^ix)] [fg'^ix) - Atgix) + l] 


-2t [fg^ix) - 3fg‘^ix) + 3tgix) + l] gfx) - ^gf) (2 -4 ^ gf) - gfx) 


= 0 . 


( 6 . 6 ) 


A. Cylindrical geometry 


In the same way, we get 


For a = 1, setting to zero the coefficient to I/x in 
dSSl), gi it) must obey the equation 


til - At + - tfg'lit) 

+ i3-8t + l2f + f)g[it) 

+2il + 3t - 3f + f) = 0. 

(6.7) 


The regular solution of this equation with the BC 51 (0) = 
— 1/8 reads 


52 (t) 


gft) 


2543 11 Sols 1 4 

2^x 32 2 x 32^25 32 ^2^x32 

^2x51 5x7 1 ^23 1 

^ 3 l-t“ 2x32(1-t)2 +^(l-t)3’ 

(6.9) 

5 X 41 X 659 72 ^ 617 2 , 17 ^3 

210 X 33 24 X 32 ^” 29 X 3^ 22 X 33 ^ 

13 X 19 4 1 ,5 _ 1 ^6 

27 X 33 23 X 32 26 X 33 

193 1 ^ 11173 1 

~ 2 X 32 1 -1 23 X 33 (1 - t)2 

4649 1 ^ 43 1 23 1 

~ 24 X 33 (1 -t)3 ^ (1 - t)4 “ ^ (1 -t)®’ 

( 6 . 10 ) 


gft) 


29 


-\—t -1 

3 12 


4 1 


etc. 


3 1 -1 


( 6 . 8 ) 


In the saturation limit A —> 00 , the requirement 
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/sat(Ka) = 2 — \/Z implies the large-Ata expansion 

79 \ 1 
na 

/17053 19765 \ 1 


/ ^ ^ AT 67 79 \ 

4atM-2-A/3+ + 


+ 


V 576 384^3/ ('«a)2 


+ - 


751423 869059 \ 1 


4608 3072^3/ (kq)^ 


+ ■ 


^sat 6 


6 


Ka 


( 6 . 11 ) 


With the aid of (13.191) . for large values of na the satura¬ 
tion value of the effective charge behaves as 


Xl^lB=H2-V3)Ka+ (-y + 


491 569 \ J_ 

6 4v^/ Ka 

43519 50329 \ 


35V3\ 

4 ) 


1 


96 


~ 0.804 Ko-f 0.655- 
0.700 „ / 1 


64v^/ («;a)2 
0.295 


+ 0 


{nay 


{nay 


O 


(na) 


( 6 . 12 ) 


The first two terms of this series have been obtained in 
Ref. [l^ . The numerical checks of the coefficients to the 
l/(Ka) and l/(Ka)^ terms are presented in Fig. [3l The 
quantity plotted in the inset is thus 


XI^Hb - 3(2 - + 


35^3^ 


(6.13) 


B. Spherical geometry 


For our purpose, the value of ^2 at t = 2 — \/3 will be 
important. Denoting 32(2 —a/S) = B, using Mathematica 
we got 

Y' 


B = -- 


IOV 3 - - -^ln2-2V31n3 

9 73 


lnl 2 - 


7 ■ 

7!. 


ln(2 - 73)-b-[ln(2 - 73)]2 

o 


ln(73 —1)-|-^ln(73-|-l) 


373 


373 


——Li2(—2 -|- Ts) -b —Li2(—6 -b 473) 

o o 

~ 1.71475. (6.19) 

In the saturation regime, the condition /sat (na) = 2 — 
73 implies the large-tea expansion 

41 \ 1 


/ 35 

tsati^a) = (2 — 73 ) + (-^ -b 


273/ 


1123 1291 


18 


A. 


1273 

sat e-“ 


- (2 - 73)-B 


1 


(Kay 

+ ... = 11^1 -• ( 6 - 20 ) 

6 Ka 

Based on (13.281) , the large-tta expansion of the saturation 
value of the effective charge is obtained in the form 

^sat 

a 


= 6(2 - 73)Ka -b (-58 -b SSTs) 


+ 


913 


1045 

273 


- 6(2 - 73)-B 


■ 1, / 

^ 1 ^ 

-b 0 1 

— 

Ka 

\Ka J 


~ l.eOStca+ 2.622- 0.089- 


o| — ) . (6.21) 

, Ka ' 


The first two terms of the expansion coincide with those 
obtained in Ref. [T^. The prefactor to the third 1 / (kq) 
term is checked against numerics in Fig. [2] 


For 0 = 2, setting to zero the coefficient to 1/x in 
dS^, gi{t) obeys the equation 

t(l — 4t + 47 — 7 ) 5 " (t) 

+(3-8t +127 + 7)7(7 

+4(1+ 3t-37+ 7) = 0. (6.14) 

The regular solution of this equation with the BC 51 (0) = 
0 reads 


51 (7 


8 4 1,81 

—I —t —y -. 

3 3 6 31 -t 


(6.15) 


Setting to zero the coefficient to l/x^ in (16.61) . we ob¬ 
tain the differential equation for g 2 (t) of type (15.151) with 
the polynomial coefficients 

P(t) = 971 - 771 - 3t - 37 - 7), (6.16) 

Q(t) = 9(1 - 7^(3 - 8t + 127 + 7), (6.17) 

i?(7 = -78 + 51t - 847 + 2527 + 19127 
-25607 + 18607 - 4307 - 3147 
+1897 - 327° + 277 (6.18) 


VII. CONCLUSION 


In this work, we have revisited the analytical results 
following from a multiple scale expansion of the non¬ 
linear Poisson-Boltzmann equation, for both cylindri¬ 
cal and spherical macro-ions. The corresponding planar 
case is analytically solvable. Three types of electrolyte 
have been addressed: symmetric ones where the co-and 
counter-ions bear the same charge in absolute value (1:1 
case), as well as asymmetric 1:2 and 2:1 situations. The 
latter two cases are not equivalent due to the non-linear 
nature of the differential equation to be solved, although 
they can yield the same Debye length. Inspecting the 
structure of the double series appearing intimates that a 
partial resummation can be performed. In doing so, and 
restricting to the 1:1 case for the sake of simplicity, the 
dimensionless electrostatic potential 4 > appears to depend 
on radial distance r through 


(j)(r) = 2 In 


( i + m \ 
Vi-/(^)y 


(7.1) 
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where x = kt, and 


Here, a is a fingerprint of geometry (more precisely, of 
curvature, with a = 0 for plates, a = 1 for cylinders, 
a = 2 for spheres) and A parameterizes the solution: 
different values of A correspond to different bare charges. 
The saturation phenomenon means that while A changes 
in some finite interval [0, Hgat], tbe bare charge varies 
between 0 and oo. More precisely, since one has (j) ^ t 
for r —>• oo, H is directly related to the effective charge of 
the macro-ion. For 1:2 and 2:1 electrolytes, relation (HID 
changes to some extent (see Eqs. (15.21) and (16.21) 1. while 
the relation between /, t and A in Eq. dia) is essentially 
unaffected. 

The planar case is such that ko —>■ oo, with f(x) = 
t = Ae~^j4:. As a consequence, our family of solutions is 
of “quasi-planar” type, which is of course quite expected 
in the limit where the macro-ion radius a is much larger 
than the Debye length 1/k. Yet, the details of this quasi¬ 
planarity are non trivial, and are such that particularly 


convenient expansion properties ensue in the asymptotic 
limit Ka ^ oo. As an illustration, we have computed 
saturated effective charges (meaning in the limit where 
the macro-ion bare charge becomes very large) where our 
scheme yields an exact expansion in inverse powers of na. 
Indeed, a careful numerical calculation of the same quan¬ 
tities from solving directly the Poisson-Boltzmann equa¬ 
tion, allows to check, term by term, the predicted expan¬ 
sion. This requires an extrapolation procedure, which 
has been presented, for extracting the saturation values 
from results that are necessarily obtained at finite al¬ 
though large bare charges. 

So far, not enough is known on the planar case for 
different asymmetries than 1:2 and 2 : 1 , so that our ap¬ 
proach cannot be generalized to such situations. What 
misses is the explicit structure of the counterpart of Eqs. 
(Tfrl) . (IQ) and (16.21) in these cases 
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